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To construct continuum stochastic growth equations for competitive nonequilibrium surface- 
growth processes of the type RD-fX that mixes random deposition (RD) with a correlated-growth 
process X, we use a simplex decomposition of the height field. A distinction between growth processes 
X that do and do not create voids in the bulk leads to the definition of the effective probability peff of 
the process X that is a measurable property of the bulk morphology and depends on the activation 
probability p of X in the competitive process RD-fX. The bulk morphology is reflected in the surface 
roughening via nonuniversal prefactors in the universal scaling of the surface width that scales in 
Peft- The equation and the resulting scaling are derived for X in either a Kardar-Parisi-Zhang or 
Edwards-Wilkinson universality class in (l-fl) dimensions, and illustrated by an example of X being 
a ballistic deposition. We obtain full data collapse on its corresponding universal scaling function 
for all p € (0; 1]. We outline the generalizations to (1 -I- n) dimensions and to many-component 
competitive growth processes. 
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I. INTRODUCTION 

Many dynamical complex physical systems in nature 
are studied by their mapping onto a suitable nonequilib¬ 
rium surface-growth problem. The dynamics of correla¬ 
tion buildup in these physical systems, and their other 
properties, can then be explored with the use of surface- 
growth methodologies. Numerous examples of such stud¬ 
ies, experimental as well as theoretical and computa¬ 
tional, come from a variety of fields such as tumor-growth 
processes [I| in cancer research, growth of cell colonies 
Q in biophysics, roughening of lipid bilayers Q in soft- 
matter biomaterials, ^namics of combustion fronts 
imbibition processes |^, film-growth processes [6|], time- 
series and market price analyses in econo-physics [7| , and 
scalabil^ and synchronization of parallel-computing sys¬ 
tem [1, Q in computer science, to give representative ex¬ 
amples. 

Large-scale properties are described within a contin¬ 
uum model by universal stochastic growth equations and 
tested with simulation models. On the theory side, the 
trouble is that simple discrete models such as SOS (solid- 
on-solid) are often not adequate to reproduce the com¬ 
plex physics of observed surface phenomena, because 
they assume only one universal process alone being re¬ 
sponsible for surface formation. Such an idealization 
does not reflect actual experimental settings where the 
observed surface phenomena may involve contributions 
from several universal processes. The continuum descrip- 
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tion of such a multicomponent growth process has not 
yet been developed. Such mixed-growth systems display 
many nontrivial properties |10l - l35l| . 

A representative example comes from an applied model 
in computer science @ when the asynchronous dynamics 
of conservative updates in a system of parallel proces¬ 
sors is modeled as a virtual-time surface that represents 
nonequilibrium processes in this system. When the load 
per processor is minimal, this dynamics belongs to the 
Kardar-Parisi-Zhang (KPZ) universality class [^. How¬ 
ever, when the load is increased to reflect real operations, 
the realistic dynamics is a competitive growth process 
that combines a universal KPZ process with a random de¬ 
position (RD) process, i.e., is of the type RD-I-KPZ (^ . 
Consequently, in order to fully understand the statistics 
of the updates and make quantitative predictions, it is 
important to know how the nonuniversal properties of 
the multi-component processes affect the universal scal¬ 
ing of a RD-f KPZ process. This is still an unsolved prob¬ 
lem of nonequilibrium surface growth science. In applied 
modeling, even if not explicitly assumed, competitive dy¬ 
namics naturally arises. Studying these systems should 
also contribute to the understanding of differences be¬ 
tween the expected and the actual scaling of rough in¬ 
terfaces, often encountered both in simulations and in 
experiments. 

By a competitive-growth process Y-l-X — alternatively 
called a two-component system or a mixed-growth pro¬ 
cess — we understand a dynamical process where process 
Y alternates with process X in accordance with the rule 
of the exclusive alternative: “either process Y (active 
with probability q) or process X (active with probabil¬ 
ity p),” is active. Here q p = 1, and Y belongs to a 
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different universality class than X. It is understood that 
in competitive surface-growth processes an event on the 
surface is triggered by only one process at a time, even if 
both mechanisms X and Y are simultaneously present. 

In this article we introduce a method by which a con¬ 
tinuum stochastic growth equation can be constructed 
for competitive growth processes. We investigate a con¬ 
nection between surface roughening and the bulk mor¬ 
phology formed during the deposition in the competitive 
growth process RD-I-X, where X is a correlated growth 
process of universal dynamics different from RD. This 
connection has been already established in simulations 
of competitive growth models and of binary growth 
of thin films [ 3 , and for diffusion-limited-aggregation 
models [13 • A new aspect of our study is to provide 
a direct theoretical link between nonuniversal proper¬ 
ties of process X, as read from the bulk, and the con¬ 
tinuum equation that underlies the observed universal 
scaling laws for the competitive RD-l-X processes. In 
this work, we derive from first principles a continuum 
equation to show that its mo del-dependent coefficients 
do reflect the bulk structure. This will lead to a distinc¬ 
tion between void-producing and simple desorption and 
adsorption processes. As discussed later, this division 
into subclasses is a necessary first step towards a theory 
of many-component processes. In particular, it explains 
variations in scale dilatations observed in RD-I-X mod¬ 
els [13, m, [m, [13, [13, d^-flll- In our analysis we use 
as an example the universal RD-I-KPZ growth process in 
(1 -|- I) dimensions and generalize our approach to other 
processes in (1 -|- n) dimensions. 

This article is organized as follows: Scaling of the in¬ 
terface width in competitive RD-I-X models is outlined 
in Sec. m where we show that the full data collapse scal¬ 
ing can be obtained in a geometric scaling given by Eqs. 
© and ®. This type of scaling was heuristically pro¬ 
posed in Ref. [131 and its explicit form was derived in 
Ref. [111. Geometric scaling confirms that the RD-I-X 
systems are in the universality class of process X [30l| . 
but such data collapse is not a dynamic finite-size scal¬ 
ing. In the remaining part of Sec. |ll] we focus on dy¬ 
namic scaling that provides a connection with stochastic 
dynamics as described by a continuum-growth equation. 
In Sec. Ell where we define the adsorption-bulk-compact 
and the dense-or-lace-bulk processes, we use a concept of 
simplectic decomposition to derive from first principles 
the stochastic growth equation for simple RD-I-X pro¬ 
cesses. Hence, we find a connection between the bulk 
morphology and the surface roughening for these pro¬ 
cesses. Results of Sec. m are discussed in Sec. IIYI where 
we demonstrate by examples that in RD-I-X processes the 
nonuniversal prefactors in Family-Vicsek universal scal¬ 
ing function are nontrivial and have connection with the 
bulk morphology. In Sec. IIVI we also give the extension 
of the approach introduced in Sec. uni to (I -I- n) dimen¬ 
sional models of two-component processes, outline a pos¬ 
sible generalization to many-component competitive pro¬ 
cesses, and discuss further developments. Conclusions 



FIG. 1: (Color on line) Typical time evolution of the inter¬ 
face width w{t) in a competitive process RD-I-X for selected 
system sizes L and activation probabilities p of process X. (a) 
Synthetic data from Eq. © when process X is in the KPZ 
universality class, (b) The collapse of the data in (a) that is 
produced by the geometric scaling {x,y) —>■ {x',y') given by 
Eqs. © and Q (shown in the figure). Here, di = 1/2 and 
d 2 = 1/3 and log(-) = logio(-). 

are summarized in Sec. El 


II. SCALING OF THE INTERFACE WIDTH 

Time-evolution of the correlation length is reflected in 
the interface width w{t) of the growing surface. Both 
the correlation length and w(t) have the same scaling 
properties. In SOS models of surfaces growing on a 
substrate of L sites, w{t) is measured as: {w'^{t)) = 
{L~^ [hk(t) — where hk{t) is the column 

height at site k at time t, and h{t) is its mean over L 
sites. The time t is measured as a number of deposited 
monolayers. The angular brackets denote configurational 
averages. For brevity of notation, we set w = y'JuF). 

In competitive growth processes that mix correlated- 
growth process X with randomness, i.e., of the type 
RD-fX, results of simulations with a flat-substrate initial 
condition at t = 0 can be summarized by the following 
two-parameter family of curves [^ : 

( ci\/t , t€[0-,ti(p)] 

w{t;L,p)=^ C 2 tP^ , t e lti{p);t 2 {L,p)) (1) 

[ , t e [t 2 iL,p);+oo], 

where p and L are parameters, and ci, C 2 , and C 3 are 
constants. The effect of the parameter p on the time 
evolution of w{t\ L,p=l) is a nonuniversal dilatation of 
time and length scales, as discussed in Ref. [131 • Because 
of these dilatations the times ti(p), when the initial RD 
transients terminate, and the times t 2 {p), when the sat¬ 
uration phases begin, have different values for different 
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curves w{t]L,p). In Eq. (IT|) (32 and a 2 are universal scal¬ 
ing exponents (the growth and the roughness exponents, 
respectively) characteristic of the universality class of 
process X. For processes RD-I-KPZ (when a = 02 = 1/2 
and (3 = 132 = 1/3) the family given by Eq. ([T]) is il¬ 
lustrated in Fig. [T^, where y = logw is plotted versus 
X = logt. In the (a:, j/)-plane, the family in Eq. ([T|) can 
be collapsed onto one curve shown in Fig. [Hd by means 
of the shift-and-scale operator G : {x,y) —>■ {x',y') [s^, 


y ^ y' 


X — Xi 1 

X 2 -Xi 132’ 
y-yi 

?/2 - ?/l ’ 


( 2 ) 

(3) 


where h = h{x,t) is the height field, x and t are the 
spatial and time coordinates, respectively, subscripts de¬ 
note partial derivatives, ri{x, t) is the white noise, and uq 
and Ao are constants. When Aq = R Eq. m is the 
Edwards-Wilkinson (EW) equation when X=EW. 
When r/Q = Aq = 0, Eq. (O defines universal RD dy¬ 
namics. 

Many simulation models of RD-I-EW and RD-I-KPZ 
growth processes [U [3(il - l^ suggest g{p) = p^ in Eq. 
(l4l) . which leads to the Family-Vicsek universal scaling 
[4g | of the average surface width w{p,t) [s^: 

= . ( 6 ) 


where x' = logP and y' = logic'. In Eqs. @ and ® 
the pairs of numbers xi = log<i(p) and yi = logw;i(p), 
and, X 2 = logt 2 (-b,p) and 1/2 = logu' 2 (-b,p), have dif¬ 
ferent values for each curve w{t;L,p). Explicitly, wi = 
w{ti{p)) and W 2 = w{t 2 {L,p)) by Eq. ([T]). Operator 
G shifts all curves in Eq. © to one position where all 
crossover points (tx, Ws&t) to saturation are mapped onto 
one point {1/(32, 3)- Subsequently, G scales the length 
(x 2 — xi)-\/r+T7^ of the correlated-growth phase in 
(x,?/)-plane for each curve to the length of 
The full data collapse obtained by G is possible because 
each curve in Eq. o has one universal footprint, where 
the initial RD transient is followed by a specific univer¬ 
sal correlation phase. Such a collapse of the data in the 
(x, ?/)-plane by geometric scaling is an illustration of the 
previously proven fact that competitive-growth pro¬ 
cesses RD-I-X are in the universality class of process X. 
It must be stressed, however, that this geometric scaling 
expressed by G does not give the universal dynamic scal¬ 
ing function that would explain the universal shape of 
the curve in Fig. [TJj in terms of finite-size scaling of the 
corresponding stochastic dynamics as described by the 
continuum model and, possibly, nonuniversal corrections 
to scaling when p ^ 1. Manipulation of Eqs. and ® 
to obtain explicitly w{t) leads back to Eq. ([T]); thus, Eqs. 
m and ([3|) do not contain any new physical information 
in addition to that already present in Eq. ©. In sum¬ 
mary to this point, the geometric scaling lacks physical 
meaning and does not connect with the Family-Vicsek 
dynamic scaling. 

A dynamic scaling hypothesis for competitive RD-I-X 
processes [s^ states that if a correlated growth X occurs 
with a constant probability p, its continuum equation 
must be invariant under the scaling 

X X , /i -)> h/g{p) , t -)> t/f{p) , (4) 

where g{p) and f{p) are arbitrary suitable functions of 
p S (0;1]. This invariance implies that f{p) — g^{p). 
When X=KPZ, the dynamic scaling hypothesis leads to 
the KPZ equation for the RD-I-KPZ mix [s^: 

ht = vof{p)hxx + {Xo/2)f^‘^{p)hl + gix, t) , ( 5 ) 


For substrates of size L, F{y) describes two limit regimes 
of evolution: F{y) ^ if p <C 1 (growth); and, 

F{y) ^ const if y 1 (saturation). In Eq. ([S]), ot and 
2 ; are the universal roughness and dynamic exponents, 
respectively, of the universality class of the correlated- 
growth process X. The scale-dilatation exponent 8 in the 
scaling prefactors in Eq. however, is nonuniversal. 
It has been observed that, in some models, d si 1 across 
universality classes and, in some other models, 0 < <5 ^ 1 
within a single universality class [13, [ 2 ^, [s^- Also, there 
are some models where the pre factors in Eq. m do not 
at all obey a power law in p |32l . |4l| . In the next section, 
we shall establish that this variation is not accidental in 
flux-conserving models but rather reflects the properties 
of the bulk of the deposited material. 


III. AB INITIO CONTINUUM EQUATION BY 
SIMPLECTIC DECOMPOSITION 

Consider aggregations where identical particles fall 
onto a substrate of L sites. On the substrate, the in¬ 
coming particles may be accepted in accordance with a 
rule that generates correlations among the sites, i.e., in 
accordance to process X. It is understood here that the 
flux of the incoming particles is uniform and time in¬ 
dependent, i.e., the average rate at which the particles 
arrive at the substrate does not vary with time and does 
not vary with position along the substrate. The corre¬ 
lated growth X occurs with probability p and competes 
with RD growth that occurs with probability q = 1 — p. 
It is understood here that the probability p remains con¬ 
stant for the entire duration of the process RD-I-X, i.e., 
the average frequency of process X does not change with 
time and does not depend on the position of the site on 
the substrate. When a particle is accepted at a site, the 
site increases its height by Ah. If, e.g., component 1 is 
RD, and component 2 is a, correlated growth in the KPZ 
universality class, their corresponding growth equations 
are 

hi,t = rii{x,t), (7) 

h2,t = i'oh2,xx + {Xo/2)hl^^ + r]2{x,t), ( 8 ) 
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where hn{x,t), n = 1,2, is the column height at x after 
time t when the component n acts alone. Assume for sim¬ 
plicity that the noise terms are of the same strength, i.e., 
rj = rji = r] 2 . In two-component growth, when both com¬ 
ponents act simultaneously together, the column height 
h{x,t) is incremented due to either of the components 
with their corresponding probabilities p and q, p + q = 1: 

Ah{x,t) = pAh 2 {x,t) + qAhi{x,t). (9) 

Here, probability p (or q) is the fraction of contributions 
to h from component 2 (or 1). For some processes this 
fraction is identical to a fraction of times when h(x) is in¬ 
cremented due to component 2 (or 1) for the times from 
0 to t. However, as explained later, this is not so for 
all processes. In Eq. (ED, Ahn is understood as “be¬ 
ing incremented due to the process n,” n = 1,2. In this 
statistical sense, Eq. ED expresses a simplectic decompo¬ 
sition of Ah(x,t) into its vertex components Ahn{x,t). 
Dividing Eq. ED by At and taking the limit At —0 
gives the equation for time rates, ht = ph 2 ,t + qhi^t, into 
which we substitute Eqs. dZD and ([5]) to obtain: 

ht — ^0ph2^xx + {Xo/2)phl^^ + {p + q)r]{x,t). (10) 

In Eq. (fTUl) . h{x,t) is the column height that rises at 
X as the result of two processes acting simultaneously 
from the beginning to time t. Here, h 2 {x,t) is the part 
of h{x,t) that was created by component 2 in this time. 
The other part was created by component 1. In other 
words, h 2 {x,t) is so far an unknown fraction of h{x,t). 
To find a relation between h and h 2 , one must consider 
nonuniversal properties of aggregation processes. 

We distinguish between the following two groups of 
surface growth processes. In one group we place all sim¬ 
ple adsorption processes with conserved flux that do not 
create voids in the bulk of the deposited material. We 
call this group adsorption-bulk-compact (ABC) growths. 
Eor example, a simple random deposition or random de¬ 
position with surface relaxation fall into the ABC cate¬ 
gory. The other group, which we call dense-or-lace-bulk 
(DOLB) growths, contains processes that are not ABC- 
type. The DOLB group includes desorption processes 
that may lead to a dense bulk as well as adsorptions that 
lead to the formation of voids. The only type of desorp¬ 
tion processes studied here are ones due to local sponta¬ 
neous desorption at the surface, not desorption processes 
where an incoming particle strikes the surface and causes 
desorption. The latter type of desorption process would 
have shadowing effects, and hence would be extremely de¬ 
pendent on the direction of the incoming particles. Note 
that the DOLB category contains flux-conserving as well 
as flux-nonconserving processes. Eor example, ballistic 
deposition and deposition to local interface minima are 
both in the DOLB group. Note that all RD universal pro¬ 
cesses are ABC-growth processes. As we show in the next 
paragraph, when component 2 is of the ABC type, proba¬ 
bilities p and q in Eq. ED express fractional contributions 
to h in terms of times, and then h 2 {x, t) = ph(x, t). This 
is not true when component ^ is a DOLB growth. 


Consider a discrete representation of events at coor¬ 
dinate X. Suppose there are t deposition events in to¬ 
tal, with ti events due to component 1, and t 2 events 
due to component 2, t = ti -\-12. In ABC growth, af¬ 
ter t events, the total column height is h = tAh, where 
contributions from components 1 and 2 are, respectively, 
hi = tiAh and h 2 = t 2 Ah. Thus, hi/h = ti/t = q and 
/ 12//1 = t 2 /t = p. Therefore, in ABC growth /12 = ph, 
and in Eqs. ED and (HUD we can identify p = p and q — q. 



FIG. 2: (Color on line) Bulk sections obtained in Monte Carlo 
simulations of surfaces, generated by: (a) random deposition 
(RD); (b) ballistic deposition (BD) that creates the bulk com¬ 
pactness c = 0.468; and, (c) the competitive RD-I-BD process 
when p = q = 1/2. Coloring indicates time intervals. In 
panel (c), within a time interval, coloring is used to differen¬ 
tiate between deposits created by RD and those created by 
BD. 
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Next consider that the component ^ is a DOLE growth 
that creates voids. Now, an individual deposition event 
due to component 2 not only increases h by A/i, but 
may also result in the creation of voids [as illustrated 
in Fig. an example of ballistic deposition]. The net 
effect is as though component 2 deposited Ah and the 
voids. Therefore, in t 2 events, its contribution to the 
column height is /i 2 = (^2 + m)Ah, where mAh reflects 
the increase in the field height due to the presence of 
voids. The component 1 is RD, i.e., ABC-type, and hi = 
tiAh. After t events, the net column height is h = hi + 
h 2 = {t + m)Ah. Thus, hi/h = ti/[t + m) < ti/t = q 
and / 12 /h = (^2 + m)/{t + m) > t 2 /t = p. The explicit 
form of these mutually complementary fractions, hn/h 
for n = 1,2, allows them to be directly measured from 
the bulk. They are, in fact, the effective probabilities 

qeS = hi/h and pes = h 2 /h (11) 

of deposition events due to components 1 and 2, respec¬ 
tively, as would result from measuring the column height. 

For some types of two-component growth with RD, 
the probability Pes can be expressed approximately as 
the power law peff = P^ HH, where the “best” expo¬ 
nent 6 can be estimated heuristically. For DOLB-type 
growth processes that produce voids, the exponent is 
i5 < 1 because PeS > p. When the component ^ is a 
DOLE growth with desorption, in the above reasoning 
one should change m —^ —to. This will give qeS > q 
and peS < p, and peS = p^ with 5 > 1. The value of 
S depends on nonuniversal particulars of the deposition 
rule of the component process 2. Therefore, in a DOLE 
growth /i 2 = Peffh, and in Eqs. m and m p = peff and 

q = 9eff- 

In general, relations h 2 {x,t) = Pesh{x,t) and p = PeS 
hold for all processes. When the correlation component is 
an ABC growth, its effective probability is identical with 
its frequency: Peff = P, provided that column-height in¬ 
crements are identical for the both processes 1 and 2. 
Thus, when process X is in the KPZ universality class, 
Eq. (flOll gives the exact stochastic dynamics for the com¬ 
petitive RD-I-X processes: 

ht = vqpIsKx + (Ao/2)p3ff/i^ -t 7?(a;, t). (12) 

When the correlation component X is a DOLE growth, 
and when the effective probability is well approximated 
by a power law the above result can be summarized 
as peff = P^, where i5 = 1 for ABC growths, and <5 1 

for DOLE growths. This result is combined with Eq. 
m to give the approximate continuum equation for the 
RD-I-KPZ mix: 

ht = VQP^^hxx + {Xo/2)p^^hl + r]{x, t). (13) 

When in Eq. ([5]) Aq = 0, the analogous reasoning gives 
the exact result for RD-I-EW dynamics: 

(14) 



FIG. 3: (Color on line) Properties of the bulk formed in com¬ 
petitive process RD+X, where X is an adsorption process that 
creates voids. When X acts alone, i.e., in the absence of RD, 
it produces a bulk with compactness c. (a) Compactness of 
the bulk c{p) as a function of the activation probability (i.e., 
frequency) p of X, plotted for selected values of c. (b) The 
effective probability Pes(p) of X for selected values of c. 


In Eq. da, we can explicitly set peff = P because all 
processes in the Edwards-Wilkinson universality class are 
ABC-type processes. When the flux particles are identi¬ 
cal, the exponent i5 = I is exact. 


IV. DISCUSSION 

Both results, Eqs. (USD, (USD and Eq. d , are in ac¬ 
cord with our former derivation that lead to Eq. ® M- 
Matching Eq. (IT^ with Eq. m gives f{p) = p^^, which 
form of /(p) was used formerly to derive the approxi¬ 
mate prefactors in Eq. ®. The inverse of the scaling 
when applied to Eqs. m and d transforms them to 
continuum equations for a “pure” correlated processes of 
p = I. Explicitly, it collapses all evolution curves w(j>,t) 
(for all L and p) either onto w{l,t) or onto a neighbor¬ 
hood of w{l,t) [^, following Eq. ®, provided the ef¬ 
fective probabilities PeS can be well approximated by the 
power-law p^. When such a ht is not possible, Eq. is 
still obeyed but then the scaling prefactors must be ex¬ 
pressed directly in terms of effective probabilities. This 
is because the factor p^ in the coefficients of Eq. d is 
only a ht to the effective probability PeS- 

Effective probabilities, dehned by Eq. d, are func¬ 
tions of the activation probability p of process X, i.e., 
the average frequency of process X in the mix RD-I-X. 
Properties of the bulk morphology created by RD-f X can 
be equivalently expressed either in terms of peff(p) or in 
terms of q^siq), because of the identities PeS + qeS = I 
and p + q = 1. Effective probability Peff of X in the mix 
RD-I-X can be expressed in various equivalent functional 
forms that may involve either the average compactness 


ht = voplshxx + vix, t). 
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c(p) [or the number density of voids of the bulk cre¬ 

ated by the mix RD-I-X or the average compactness c [or 
the number density v of voids] of the bulk created by the 
process X acting alone (i.e., in the absence of RD when 
P = 1): 


Pes{p) = ^ = i’{p)+pc{p) = 1 - qc{p), (15) 

h m + t 

where c{p) = t/{m + t), v(p) = m/{m + t), and c(p) + 
v{p) = 1; 


Pes{p) 


h2 _ m + t2 _ p _ p 
h m + t c + pv p + qc’ 


(16) 


where c = < 2/^21 v = m/h 2 , p = ^ 2 /^, and c + v = 1. 
Combining Eq. (ITKll with Eq. (I16|) gives the compactness 
c(p) of the RD-I-X bulk as a function of the activation 
probability p of process X: 


c{p) = 


c 

c-\- pv 


(17) 


in Refs. [ 2 ^, |3T| for RD-I-BD models, does not produce 
data collapse at all. The RD-I-BD model when BD is the 
next-nearest-neighbor (NNN) sticking rule [l^ provides 
an example where pes{p), and thus the nonuniversal pref¬ 
actors f{p) and g{p) in Family-Vicsek universal scaling, 
cannot be expressed by a power law p^ [s^ . In this sys¬ 
tem the surface roughening obeys power laws in effective 
probability that incorporates either the compactness or 
the void density of the bulk, resulting in excellent data 
collapse of w{p,t), similar to that seen in Fig. mmi. 

The approach introduced here by the example of a 
KPZ processes, can be applied to a broad range of 
stochastic growth models RD-f X, where component 2 can 
be any isotropic growth in (1 -|- n) dimensions: 

h2,t(.x,t)=D{h2) + m(.x,t), (19) 

where x is n dimensional, and the operator D represents 
only local interactions [^. In the general case, Eq. m 
is written as ht = Peff^ 2 ,t + and combined with 

Eqs. © and (USD to find for the competitive growth 


Equations (HID, (USD, and (ED show that c{p) and Pes{p) 
can be easily measured in experiment as well as in simu¬ 
lations: All that is needed is to measure the average den¬ 
sity of voids in a sample cross section of the bulk when 
the correlated process X acts alone (i.e., in the absence 
of RD whenp =1). Equations (fTHll and (flTll . plotted for 
several values of c in Fig. |3l show that neither the com¬ 
pactness c(p) nor the effective probability Pesip) follow 
a power law in p when c < 1. 

The dynamic scaling hypothesis for RD-I-X processes, 
Eq. O, can be reinstated by expressing f{p) explicitly 
in terms of Peff: f{Pes) = g'^iPes), since Eq. (|T6l) can 
be inverted to give p{pes). Repeating the steps outlined 
in Ref. [s^ gives the following generalization of Family- 
Vicsek scaling for the surface roughness: 

L^a / ^ \ 

w'^ip,t) = ^^ Trd+X yf{PeS)j^J , (18) 

where Erd-i-x(’) describes the three regimes of the evo¬ 
lution seen in Fig. [Do. The effect of the nonuniversal 
prefactors f{p) in Eq. (IT^ is a dilatation of length and 
time scales, as discussed in Ref. [s^. The physical mean¬ 
ing of p is that of a noise-tuning parameter. 

In Fig. |4]we give an example of the exact scaling where 
nonuniversal prefactors in Eq. m are directly expressed 
by PeS via the substitution p^ —>■ pfff(p) = \/ f{p) for 
the RD-fBD model when ballistic deposition (BD) is the 
nearest-neighbor (NN) sticking rule [^. Here, the effec¬ 
tive probability depends on both p and the mean com¬ 
pactness c{p) of the bulk formed in the RD-I-BD process, 
given by Eq. ED ( 41 I I . The excellent data collapse in the 
full range oi p G (0; 1], seen in Fig. 21 can be contrasted 
with Fig. 5 of Ref. [30| that shows only an approximate 
data collapse for the same system with the best fit expo¬ 
nent 6 ~ 0.41 in Eq. ([HD. It needs to be said explicitly 
that the scaling where <5 = 1/2 in Eq. ([HD, proposed 


ht{x, t) = p^sb{pesh) + riix, t) , (20) 

where 77 = (1 — Peff)? 7 i + PeS'il 2 , and the noise strengths 
may be different. Equations ED and ( 1 ^ represent the 
same universality class since the multiplication by peff 
does not modify local interactions: Peff affects the noise 
strength and the gradient of the height field, but does 
not generate any new terms other than those already 
given by operator D. Hence, if a correlated growth be¬ 
longs to a given universality class, its mix with RD will 
remain in the same universality class. Elementary calcu¬ 
lations show that Eq. (EOD is invariant under the scaling 
g{p)h{x,t) = h'{x,t' = f{p)t). If g{p) = Pes{p) and 
f{p) = Peff(p)i if the noise strengths are the same, 
this scaling maps the universal dynamics (1201) of RD-I-X 
onto the universal dynamics of X. In this case the in¬ 
variance implies g(j))w{p,t) = w'{f{p)t), where w'{-) has 
universal scaling properties of the process X. When X is 



FIG. 4: (Color on line) Scaled time-evolution nP{p,t) in the 
RD-I-BD model. In this example, the scaling function f(p) = 
g^(p) explicitly incorporates the compactness c(p) of the bulk 
formed in the RD-I-BD process. Here, L = 500, 2a = 1, and 
averaging was performed over 400 surface configurations, i.e., 
independent simulations. 
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either in the KPZ or in the EW universality class, and if 
additionally Peff ~ ^ we recover Eq. 

When both the RD and the correlation component 2 
have deposits of unit height, when peff ~ we have 
(5 = 1 if component 2 is of the ABC-type; and, d ^ 1 if 
it is of the DOLE type. In the latter case, the value 
of the exponent <5 is specific to component 2. When 
Pes incorporates explicitly bulk properties, the scaling is 
g{p) = pfjj(p), where the new scale-dilatation exponent (5 
is obtained from the slope of lnr(;^(p) plotted vs Inpeff(p) 
at saturation. In DOLE growth with voids, PeS can be 
determined by measuring the mean density of voids in 
the bulk (Figs. [2]and|4l). Similarly, in DOLE growth 
with desorption, peff is connected to the mean fraction of 
the removed material (or flux) [4l|. 

The analysis presented here explains scaling results of 
the followiM mix ed-g rowth models in (1 -|- 1 ) dimensions. 
Model A [ 13 , HH, l30| : component 2 is RD with surface 
relaxation. Model B : component 2 simulates a depo¬ 
sition of a sticky non gran ular material of variable droplet 
size. Model component ^ is the NN sticking 

rule of ED. Model D |22l : component .Sis a deposition 

of Poisson-random numbers to the local surface minima. 
Models A and B are ABC growths in the EW universal¬ 
ity class, where Peff = P and 5 = 1 (Fig. [5]). Models C 
and D belong to the KPZ universality class. Model C is 
an example of DOLE growth with voids, with a 53.2% 
void density in the bulk when p = 1 , and in this case 
5 ~ 0.41 < 1. Model D is a DOLB-type growth that 
produces a compact bulk but component 2 is flux non¬ 
conserving, and here <5 Ri 1. Extensions of Models A and 
C to (1-1- n) dimensions [2l|, n = 2,3, yield results that 
conform to our theoretical predictions of Peff ~ P^ with 
(5 7 ^ 1 for mixing RD with DOLE processes, and d = 1 
for mixing RD with ABC processes. Additional exam¬ 
ples include cases [ 2 ^ when component ^ is a restricted 
Kim-Kosterlitz solid-on-solid model (where RD-I-X is 
in KPZ universality class), and when it simulates a con¬ 
served restricted SOS growth of Kim et al. [i^. In the 
latter case the process RD-I-X is in the Villain-Lai-Das 
Sarma universality class [l^, IH, . 

An interesting lattice simulation model has been re¬ 
cently considered by Banerjee et al. [35l| in an attempt 
to describe a realistic sedimentation. The Banerjee et 
al. model is a competitive growth process that has three 
component processes: one RD process and two DOLE 
processes, where one DOLE process is ED with the NN 
sticking rule and the other DOLE process is ED with the 
NNN sticking rule. In the language of our study, the over¬ 
all process is the RD-I-X process, where X=X1-|-X2 and, 
as the convex linear combination of two KPZ processes, 
X is in KPZ universality class. Accordingly, this system 
should obey the scaling law of Eq. m in the effective 
probability of the combined process X. Time-evolution 
plots of the surface roughness in Ref. suggest such 
scaling. 

The extension of the approach presented here to other 
competitive growth processes may provide a tool to un- 
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FIG. 5: Scaled widths at saturation w vs the parameter 1/p'*. 
Panels (a) and (b) are for models A and B, respectively. Pan¬ 
els (c) and (d) are for models C and D, respectively. Reference 
lines have slope 1. Data are scaled with the a values shown. 


derstand the observed dynamics of surface growth. Re¬ 
alistic systems may involve many component processes, 
some of which may be dominant. Within our formalism 
a departure point may be a generalization of Eq. ®: 

A/i(f , t)=^ {x, t), (21) 

k 

where the summation is over all contributing processes, 
and is the column-height increment due to the fcth 

process. In a first approximation, component processes 
are not explicitly correlated. Each process is encoun¬ 
tered with the activation probability or frequency pk, 
— 1) contributes to the growth with an ef¬ 
fective probability p^^g , J2kPeS = 1- In tbe trivial case 
when all components are ABC-type models with the unit 

(k) 

mean deposit height we have p^^ = pk. For a DOLB-type 

(k) 

growth the p/g will have to be determined. Depending 
on the model, p/g can be estimated by analyzing the 
growth when process k acts alone, and measuring either 
the mean bulk density or the mean fraction of the de¬ 
tached material or both [4l|. Simplectic decompositions 
like the one proposed in Eq. m have a long history of 
applications in many diverse fields and are the precursors 
of probability measures. 

Stochastic theory of multicomponent competitive far- 
from-equilibrium surface-growth processes is a newly 
emerging topic in statistical physics. During the re¬ 
cent two decades, in addition to model-specific simulation 
studies of two-component (either RD-I-EW or RD-I-KPZ 
lattice) growth models, a special case of RD-fRD has 
been considered both in a theoretical mean-field ap¬ 
proach and in simulations [d^. The absence of a con¬ 
sistent continuum theory for the RD-I-X mix hindered 




































scale-invariance studies of more complex systems such 
as EW-I-EW or KPZ-I-KPZ or EW-I-KPZ, or more gen¬ 
eral three-component systems (such as, e.g., those in Ref. 
[s^). Understanding the dynamics of real physical sys¬ 
tems calls for more realistic models that would go be¬ 
yond a one-process theory of kinetic roughening. For 
example, in realistic modeling of ion-bombardment ex¬ 
periment, where shadowing effects matter, the mecha¬ 
nism of ion desorption must be (at least) accompanied 
by ion diffusion along the substrate as well as by random 
deposition. The construction of a continuum stochastic 
growth equation, outlined in this article, will be helpful 
for future studies of two- and three-component competi¬ 
tive nonequilibrium growth systems. 

V. CONCLUSION 

In summary, we derived continuum stochastic-growth 
equations and the resulting scaling for competitive 
RD-fX growth processes. The RD-I-X growth processes 
show that mo del-dependent prefactors in universal scal¬ 
ing laws can be linked with the bulk morphology and 
determined from bulk structures. This necessitates the 


distinction between the adsorption-bulk-compact (ABC) 
and the dense-or-lace-bulk (DOLE) growth processes in 
dynamic-scaling analysis of competitive mixed-growth 
models. For competitive systems, the activation prob¬ 
ability of process X, i.e., its frequency in the RD-fX mix, 
alone does not provide sufficient information to correctly 
describe their dynamics. The essential physical nonuni- 
versal parameter here is the effective probability, p^ff, of 
X. The bulk morphology allows one to obtain peff for 
either experimental or computational studies. Further¬ 
more, Peff is parametrized by the activation probability. 
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